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ABSTRACT 

We use cosmological hydrodynamic simulations with unprecedented resolution to study the forma- 
tion of primordial stars in an ionized gas at high redshifts. Our approach includes all the relevant 
atomic and molecular physics to follow the thermal evolution of a prestellar gas cloud to very high 
densities of ~ 10^*cm~^. We locate a star- forming gas cloud within a reionized region in our cos- 
mological simulation. The gas cloud cools by HD line cooling down to a few tens Kelvin, which is 
lower than possible by H2 cooling only. Owing to the low temperature, the first run-away collapse is 
triggered when the gas cloud's mass is ^ AOMq. We show that the cloud core remains stable against 
chemo-thermal instability and also against gravitational deformation throughout its evolution. Con- 
sequently, a single proto-stellar seed is formed, which accretes the surrounding hot gas at the rate 
M ^ lO^^M0/yr. We carry out proto-stellar evolution calculations using the inferred accretion rate. 
The resulting mass of the star when it reaches the zero-age main sequence is Mzams ~ 40Mq. We 
argue that, since the obtained Mzams is as large as the mass of the collapsing parent cloud, the 
final stellar mass should be close to this value. Such massive, rather than exceptionally massive, 
primordial stars are expected to cause early chemical enrichment of the Universe by exploding as 
black hole- forming super/hypernovae, and may also be progenitors of high redshift 7-ray bursts. The 
elemental abundance patterns of recently discovered hyper metal-poor stars suggest that they might 
have been born from the interstellar medium that was metal-enriched by supernovae of these massive 
primordial stars. 

Subject headings: stars:formation 
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1. INTRODUCTION 

Recent observations of high redshift galaxies 
(Kashikawa et al. 2006), quasars (Fan et al. 2006), and 
gamma-ray bursts (Totani et al. 2006) provide rich in- 
formation on the ionization structure of the intergalactic 
medium (IGM) at z ^ 6 — 7. The observations indicate 
that a large fraction of the IGM was already ionized 
when the Universe was less than one billion years old, 
but suggest also that a substantial fraction was still 
neutral. Combined with polarization measurements of 
the cosmic microwave background (CMB) radiation 
by the WMAP satellite (Spergel et al. 2007; Page et 
al. 2007), it is inferred that cosmic reionization was 
extended and lasted for several hundred million years 
(see, e.g. McQuinn et al. 2006). 

The standard cosmological model based on cold dark 
matter and dark energy predicts that the first stars form 
in low-mass halos at high redshifts (Couchman & Rees 
1986; Yoshida et al. 2003). Reionization subsequently 
proceeds in a patchy manner, whereby Hii regions de- 
velop around individual sources and eventually percolate 
(Gnedin & Ostriker 1997; Miralda-Escude et al. 2000; 
Sokasian et al. 2003, 2004; Furlanetto et al. 2004, 2006; 
Zahn et al. 2007). Since the sources, massive stars, are 
expected to emit ionizing photons for only a few to ten 
million years, the reionized IGM will start recombining 
and cooling when the stars die. We study the formation 
of primordial stars in such reionized regions. 
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While there have been a number of numerical stud- 
ies of the formation of primordial stars in a neutral gas 
(Omukai & Nishi 1998; Abel et al. 2002; Bromm et al. 
2002; Yoshida et al. 2006, hereafter Y06), the evolu- 
tion of a primordial gas cloud during/after reionization 
has not been examined in detail. It is well-known that 
the thermal evolution of a recombining gas differs sig- 
nificantly from initially neutral primordial gas (e.g., Na- 
gakura & Omukai 2005; Johnson & Bromm 2006). More- 
over, it has been suggested that cooling by HD molecules 
leads to the formation of low-mass, possibly near solar- 
mass, primordial stars (Uehara & Inutsuka 2000; Naka- 
mura & Umemura 2002). It is important to determine 
whether or not such small mass gas clumps are formed 
in a reionized IGM at high redshifts. 

Earlier (Yoshida et al. 2007; hereafter Y07), we used 
a large cosmological simulation to study the evolution of 
early relic Hii regions until second-generation gas clouds 
are form. In this Letter, we further explore the evolution 
of these prestellar gas clouds. We investigate in detail the 
structure of such gas clouds and examine whether they 
become chemo-thermally unstable to fragmentation into 
smaller clumps. We then compute the gas mass accretion 
rate and use it as an input to a proto-stellar calculation. 

2. PHYSICAL CONDITIONS AROUND THE PROTOSTAR 

We use the simulation outputs of Y07. Briefly, they 
carried out radiation-hydrodynamics calculations of the 
formation and evolution of early Hii /Hem regions 
around the first stars. They showed that the ionized 
gas cools and re-collapses in a growing dark matter halo 
within a hundred million years, to form so-called second 
generation star(s). They also examined the impact of an 
external far-UV background radiation field. 
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Fig. 1. — Top: Radial temperature profile as a function of en- 
closed gas mass. We explain the characteristic features as follows; 
(A) Efficient molecular hydrogen production in the recombining gas 
enables rapid cooling to a few hundred Kelvin; (B) HD molecules 
are formed and the gas cools further below 100 K; (C) heating by 
the CMB floors the gas temperature, and the cloud starts to con- 
tract gravitationally; (D) three-body reactions kick in and the core 
becomes fully moleculax; (E) the core is already optically thick to 
H2 lines but cools by collision-induced emission; (F) the core be- 
comes optically thick to continuum. Bottom: Instantaneous gas 
mass accretion rate as a function of enclosed gas mass. The solid 
circle indicates where the enclosed gas mass exceeds the locally es- 
timated Bonnor-Ebert mass. The dot-dashed line is the accretion 
rate around the first generation star in the simulation of Y06. 

We use the output of the fiducial run of Y07, in which 
a second generation star forms at 2; = 16 in a halo with 
mass ~ 10 ''Mq, whose virial temperature is ~ 10^ K. In 
addition to the atomic and molecular physics described 
in Y07, wc have also implemented comptitation of local 
optical depth to continuum radiation using the Planck 
opacity table of Lenzuni et al. (1991). The contin- 
uum opacity is used to evaluate the net cooling rate by 
collision- induced emission. 

Figure 1 shows the radial profiles of gas temperature 
and gas mass accretion rate in the star-forming gas cloud. 
We use the final output of the simulation when the cen- 
tral density is lO^^cm"^. The density profile at this time 
is close to a single power-law of p oc r~^'^ over nearly 
20 decades in density. The characteristic features of the 
temperature profile can be understood by appealing to 
various atomic and molecular processes, as denoted in the 
figure and explained in the caption. Radiative cooling of 
a molecular gas is halted by effective heating from CMB 
photons that have a temperature of Tcmb = 2.728(1-1-2;) 
at redshift z. The minimum temperature of the molec- 
ular gas cloud is limited to this value. We calculate the 



Bonnor-Ebert mass using the density, temperature and 
molecular fraction profiles; the cloud is dynamically un- 
stable at a mass scale of ~ 4OM0. This is the mass for the 
first run-away collapse of the star-forming cloud. Note 
that it is essentially the Jeans mass at the temperature 
minimum. 

The gas mass accretion rate shown in the bottom panel 
of Figure 1 reflects the temperature structure. Note that, 
overall, the accretion rate is much smaller than that for 
the first star simulated by Y06 (dot-dashed line). This 
can be understood by the fact that the mass accretion 
rate is roughly M oc /G, where Cg is the sound speed. 
HD cooling enables the gas cloud to cool to lower tem- 
peratures than is possible by H2 cooling, and hence the 
accretion rate within the cloud is smaller. 

Wc do not consider cooling by heavy elements in our 
simulation, assuming that the gas is not chemically en- 
riched. However, the overall thermal evolution will not 
differ much from that shown in Figure 1, unless a signifi- 
cant amount of dust is present, because cooling by metal 
lines is efficient only at low densities. HD cooling alone 
can bring the gas to Tcmb , and the gas temperature can- 
not fall below this even with atomic metal cooling. 
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Fig. 2. — Growth parameter Q = tdyn/tg versus gas density. We 
use an output when the central density is lO^^cm"^ and show the 
distribution of the gas elements in this plane as grey-contours. All 
of the gas elements are in the stable region Q < 1, where the dy- 
namical time is shorter that the growth time of the perturbations. 



3. THERMAL INSTABILITY AND FRAGMENTATION 

An important question is whether or not the gas cloud 

becomes thermally unstable owing to efficient molecule 
production and cooling which operates at densities n > 
lO^cm^"*. We calculate the growth rate of isobaric per- 
turbations following Omukai & Yoshii (2003). We solve 
the dispersion relation for perturbations in temperature, 
density, and molecular fraction, {5T,Sp,Sf) oc exp(wi), 
under isobaric conditions. We define the growth param- 
eter as the ratio of the characteristic growth time scale 
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tg of the perturbations to the dynamical time tdyn: 

^ , ^dyn 

— Wtdyn — — • 

tgr 



(1) 



Figure 2 shows Q as a function of local gas density. 

Clearly, all the gas elements are in the stable region 
Q < 1, where the dynamical time is shorter than the 
growth time scale. We thus conclude that the chemo- 
thermal instability does not operate in this regime. Note 
that, because of the harder effective equation of state 
than the initially neutral primordial gas in the relevant 
density and temperature range, gravitational deforma- 
tion is loss effective for the gas cloud. Wc followed 
Omukai ct al. (2005) and checked that the cloud is also 
stable against gravitational deformation. Our simulation 
indeed shows that the cloud core remains stable against 
fragmentation, as expected. 

4. PROTO-STELLAR EVOLUTION 

We carry out proto-stellar evolution calculations to 
study the growth of the protostar in detail. We employ 
the scheme described by Stabler, Shu, & Taam (1980) as 
modified by Stabler et al. (1986) and Omukai & Palla 
(2003). We use the mass accretion rate shown in Figure 1 
as input to the protostellar evolution code, assuming that 
the gas is accreted in an approximately spherical manner. 
As is shown in Figure 1 , the central part is accreting the 
surrounding gas at a rate ~ lO~^M0/yr and thus a star 
with mass ~ lOM© will form within 10^ years. However, 
the final stellar mass is determined by processes such 
as radiative feedback from the protostar. Wc treat the 
evolution of a protostar as a sequence of a growing hy- 
drostatic core with an accreting envelope. The ordinary 
stellar structure equations are applied to the hydrostatic 
core. The structure of the accreting envelope is calcu- 
lated under the assumption that the flow is steady for a 
given mass accretion rate. 

Figure 3 shows the resulting evolution of the proto- 
star. After a transient phase and an adiabatic growth 
phase at M* < lOM©, the protostar enters the Kelvin- 
Hclmholtz phase and contracts by radiating its thermal 
energy. When the central temperature reaches 10*K, hy- 
drogen burning by the CNO cycle begins with a slight 
amount of carbon synthesized by helium burning. This 
phase is marked by a solid circle in the figure. The energy 
generation by hydrogen burning halts contraction when 
the mass is 35Mq and its radius is ~ 2.8 solar radii. 
Soon after, the star reaches the zero-age main sequence 
(ZAMS). The protostar relaxes to a ZAMS star within 
about 10^ years from the birth of the protostellar seed. 
Accretion is not halted by radiation from the protostar 
to the end of our calculation. 

It is important to point out that the mass of the par- 
ent cloud from which the star formed is Mdoud ~ 40Mq. 
The final stellar mass is likely limited by the mass of 
the gravitationally unstable parent cloud. We thus ar- 
gue that primordial stars formed from an ionized gas 
are massive, with a characteristic mass of several tens 
of solar masses, allowing overall uncertainties in the ac- 
cretion physics and also the dependence of the minimum 
gas temperature on redshift (see next section) . They are 
smaller than the first stars formed from a neutral gas, but 
are not low-mass objects as suggested by earlier studies. 
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Fig. 3. — Evolution of the proto-stcUar radius and the mass (solid 
line) . The solid circle marks the time when efficient hydrogen burn- 
ing begins. The dotted line shows the mass (and radius) growth 
which is calculated under the assumption that a larger amount of 
gas than the parent cloud can be accreted. For reference we also 
show the result from Y06 for a first generation star that forms in 
a neutral gas cloud. 

5. DISCUSSION 

Recent theoretical studies of the formation of the first 
stars indicate that they were likely rather massive, with a 
characteristic mass greater than a hundred solar masses 
(Abel et al. 2002; Omukai & Palla 2003; Y06; Gao et 
al. 2007). We have explicitly shown, for the first time, 
that primordial stars formed from a reionized gas have a 
smaller mass (r-^ AOMq) than these first stars. The low 
temperature of the gas cloud owes to HD cooling, and 
the correspondingly small collapse mass and low accre- 
tion rate lead to the smaller final stellar mass. We note 
that various proto-stellar feedback effects such as mass 
outflow, which wc do not take into account, could briefly 
shut off accretion (Machida ct al. 2006), and then the 
final mass would be even lower. 

Intriguingly, the elemental abundance patterns of hy- 
per metal-poor stars discovered by Christlieb et al. 
(2002) and Frebel et al. (2005) indicate that early metal- 
enrichment was caused by supernova explosions with a 
progenitor mass of ^ 25 Mq (Umeda & Nomoto 2003; 
Iwamoto et al. 2005). Massive primordial stars formed 
in reionized regions, such as those studied here, may be 
responsible for supernova explosions in the early Universe 
(Tumhnson et al. 2004; Greif & Bromm 2006). 

We have shown that the minimum temperature of the 
gas cloud is limited by the CMB, and thus the thermal 
evolution is slightly different if the gas cloud collapses 
late at a lower redshift. We ran simulations varying 
the collapse redshift (in practice, the CMB temperature) 
to examine this effect. The simtilations show that the 
mass scale for the first runaway collapse becomes smaller, 
down to Mcioud ~ IOM0, owing to the decrease in the 
minimum gas temperature. While this might seem to 
suggest a lower characteristic mass for primordial stars 
formed at lower redshifts, the effect of photo-dissociating 
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radiation in the far-ultra violet (FUV) needs to be con- 
sidered carefully in this context, because a very weak 
FUV radiation field is enough to destroy the H2 and HD 
molecules in a primordial gas cloud. 

We carried out one-zone calculations of a collapsing 
primordial gas as in Nagakura & Omukai (2005), but 
with photodissociation of H2, HD, and H~. The temper- 
ature evolution for primordial clouds with different FUV 
strength is shown in Figure 4. We identify a critical flux 
of Go — 10"^ in terms of the Habing parameter (Habing 
1968), or equivalently 3 x 10~^^erg s~-^cm~^Hz^^str~^, 
above which the formation of (and hence cooling by) H2 
molecules is limited by photo-dissociation, and the gas 
cloud starts contracting gravitationally before its tem- 
perature is lowered to < 100 K by HD cooling. The 
critical intensity is likely to be reached before the com- 
pletion of global reionization (e.g., Haiman et al. 1997; 
Yoshida et al. 2003). It is also expected that FUV 
radiation from a nearby source exceeds this limit. For 
example, a massive star with a few tens solar-masses 
emits photons in the Lyman- Werner bands at a rate of 
-^LW ~ lO^^erg s~^Hz~^. For a typical physical condi- 
tion in an early relic Hii region (see Y07), HD forma- 
tion will be suppressed if another massive star is located 
within a radius of rdiss. — 80 (riHa/lO^^cm"^)"'^/-'^^ pc 
even if we account for gas self-shielding. The number 
of stars formed by HD cooling is likely limited within a 
relic Hii region, and thus the formation efficiency of these 
stars appears to be low. 

It is also worth discussing the effect of metallicity. 
Atomic metal line cooling at low densities will not signif- 
icantly affect the evolution of the collapsing gas clouds, 
because Tcmb determines the minimum gas temperature 
even in the presence of metal cooling. On the other hand, 
the existence of dust may play an important role in signif- 
icantly changing the characteristic mass of collapsing gas 
clouds. Schneider et al. (2006) show that thermal dust 
emission in such clouds can result in low temperatures 
at high gas densities, lowering the collapse mass scale to 
below one solar mass. Therefore our conclusion remains 
robust unless a significant amount of dust is present in 
the gas clouds. 

Finally, there is an interesting possibility of observ- 



ing the death of massive primordial stars as gamma-ray 
bursts. Tominaga et al. (2007) recently suggested that 
primordial stars with mass ~ 40Mq trigger long-duration 
gamma-ray bursts and faint supernovae, such as those 
observed in the nearby Universe (e.g., Fynbo et al. 2006). 
Future observations of bright gamma-ray bursts and the 
absorption lines in their radiation spectra will provide in- 
valuable information on the interstellar medium at high 
redshifts (Inoue, Omukai, & Ciardi 2005). 
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Fig. 4. — Temperature evolution of an initially ionized primordial 
gas with FUV radiation. The FUV radiation flux is parametrized 
by the Habing parameter Go, which corresponds to the mean in- 
tensity of ~ 3 X 10~^''Goerg s~lcm~^str~^Hz~^ in the LW wave- 
length range. The initial number density, temperature and ioniza- 
tion degree are lcm~^, 7000K, and 0.1, respectively. 
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